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Hard  limiting  before  pulse  compression  or  correlation  processing  is  a 
common  approach  to  the  CFAR  (constant  false  alarm  rate)  problem,  and  it 
offers  a  good  and  simple  solution  in  a  single -target  or  scarce-target  envi¬ 
ronment.  With  the  advent  of  radars  with  a  large  time-bandwidth  product  the 
possibility  arises  that  expanded  radar  returns  due  to  multiple  targets  of  in¬ 
terest  may  overlap  very  largely  or  entirely  but  still  may  be  sufficiently  sep¬ 
arated  to  be  resolved  after  receiver  processing.  In  this  case  the  compressed 
pulses  cannot  attain  full  amplitude  at  the  processor  output  even  if  the  signal- 
to-noise  ratio  at  the  input  is  very  high;  tills  phenomenon  is  known  as  capture 
and  small  signal  suppression.  The  purpose  of  this  report  is  to  exhibit  that, 
in  addition  to  compressed  target  responses  of  reduced  magnitudes,  false 
targets  may  be  generated  with  apparent  amplitudes  of  the  eame  order  or  ex¬ 
ceeding  those  of  legitimate  targets.  Spurious  target  generation  in  the  case  of 
chirp  radar  has  been  known  for  some  time.  The  theory  has  been  extended  to 
maximum-length  linear  shift-register  codes  which  are  used  as  modulation 
functions  of  pulse-compression  and  phase -coded  CW  radars.  It  is  found  that 
a  single  pair  of  radar  returns  coded  in  this  manner  is  subject  to  capture  only 
and  not  to  false  target  generation.  Surprisingly,  however,  the  addition  of  a 
third  expanded  signal  produces  a  spurious  response.  This  generation  of  a 
false  target  should  be  taken  in  account  when  the  dynamic  range  of  future 
phase-coded  radars  using  linear  shift -register  codes  is  specified,  in  particu¬ 
lar  if  the  radar  is  designed  for  automatic  track  and  raid-size  determination. 

General  formulas  were  derived  to  predict  the  effects  of  capture  and 
false  target  generation  as  a  function  of  the  signal  energy  distribution  and 
relative  phasing  before  entering  the  limiting  device.  The  formulas  were 
evaluated  numerically,  with  the  results  being  presented  in  the  form  of 
computer-generated  plots. 
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CAPTURE  AND  SPURIOUS  TARGET  GENERATION  DUE  TO  HARD 
LIMITING  IN  LARGE  TIME -BANDWIDTH  PRODUCT  RADARS 


SUMMARY 

This  report  contains  the  results  of  a  study  of  the  effects  of  limiting  combined  with 
various  types  of  pulse  coding  in  the  suppression  or  "capture"  of  real  targets  and  the 
generation  of  false  targets.  In  the  study,  the  radar  return  is  described  by  a  complex 
signal  vector  which  is  modulated  in  amplitude  and  in  phase.  The  limiter  generates  an 
output  vector  in  phase  with  the  input  vector,  but  with  a  constant  (unit)  amplitude. 
Throughout  this  report  is  is  assumed  that  the  limiter  input  waveform  is  the  sum  of  two 
or  three  phase-coded  signals  of  the  same  kind  but  with  different  delays.  The  component 
signals  are  assumed  to  overlap  entirely.  Their  amplitudes  and  the  rf  phases  of  their 
carriers  may  be  arbitrary.  It  is  assumed  that  the  signal-to-noise  ratio  is  high  and  that 
the  beat  products  coming  out  of  the  limiter  can  therefore  be  predicted.  It  is  shown  in  the 
different  sections  of  the  report  how  the  beat  products  may  interfere  with  the  legitimate 
target  returns,  thus  causing  an  apparent  amplitude  change  (capture  effect)  and  how  they 
may  combine  and  form  new  signals  of  the  same  kind  as  radiated  by  the  radar,  thus  caus¬ 
ing  a  false  target  response  (spurious  target,  ghost  target). 

In  the  first  two  sections  of  the  report  the  model  and  the  assumptions  are  explained 
and  are  related  to  radar  designs.  Examples  show  how  limiting  takes  place  in  radar  re¬ 
ceivers.  There  are  two  cases:  intentional  limiting,  to  obtain  CFAR  (constant  false 
alarm  rate)  or  to  reduce  equipment  complexity,  and  accidental  limiting,  which  occurs  if 
the  radar  receiver  is  overdriven  by  large  clutter  returns  or  electronic  interference. 

The  sections  following  the  first  two  sections  are  devoted  to  the  analysis  of  the  lim¬ 
iter  output  if  the  input  consists  of  two  or  three  mutually  delayed  expanded  radar  signals. 
The  sections  are  the  following: 

1.  General  Theory  of  the  Hard  Limiting  of  Two  Input  Signals.  In  this  section  general 
formulas  are  derived  and  discussed.  The  nonlinear  relationship  between  limiter  output 
and  the  instantaneous  phase  difference  of  the  input  signals  is  developed  into  a  Fourier 
series.  The  Fourier  coefficients  c„  are  calculated  for  the  order  />  in  the  range  between 
-14  and  +15.  The  results  are  plotted  and  printed  out  for  a  number  of  parameter  choices. 

2.  Limiting  of  a  Pair  of  Constant -Frequency  CW  Sinusoids.  This  section  may  be  of 
interest  to  the  designers  of  CW,  pulsed  CW,  or  pulse  doppler  radars.  The  theory  com¬ 
pares  favorably  with  a  bench  test. 

3.  Limiting  of  a  Pair  of  Linearly-Frequency-Modulated  Signals.  It  is  shown  that  in 
addition  to  small  signal  suppression  there  is  a  false  target  generation  effect.  After  pulse 
compression  an  array  of  false  targets  appears  to  both  sides  of  the  true  target  returns. 

4.  Limiting  of  a  Pair  of  Zero-Pi -Phase -Modulated  Signals.  All  components  of  the 
limiter  output  may  be  identified  with  Images  of  the  original  input  signals.  The  smaller 
signal  will  be  captured  to  an  amount  depending  on  the  intensity  ratio  before  limiting  and 
on  the  carrier  rf  phase.  The  capture  effect  is  minimised  if  the  carriers  are  90  degrees 
out  of  phase.  There  is  no  evidence  of  any  false  targets. 
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5.  General  Theory  of  the  Hard  Limiting  of  Three  Input  Signals.  The  two-signal 
theory  is  extended  to  the  three-signal  case  in  a  straightforward  manner.  A  two- 
dimensional  Fourier  series  is  used  to  express  the  nonlinear  relation  between  output  and 
input  quantities. 

6.  Limiting  of  a  Triplet  of  Linear  Coded  Sequences.  It  is  shown  that  the  limiter 
output  is  composed  in  this  case  of  four  coherent  signals.  Three  signals  are  identical 
with  the  input  signals,  and  they  are  the  true  target  responses;  but  the  fourth  signal  has  a 
pseudo-random  delay,  and  it  is  a  false  target  or  ghost  target.  The  theory  is  confirmed 
by  the  results  of  a  computer  simulation.  (To  be  exact  regarding  the  history  of  this  study, 
the  computer  simulation  was  made  first,  and  the  theoretical  explanation  for  the  false 
target  generation  effect  was  sought  and  found  afterward.)  Pseudo-three-dimensional 
plots  show  the  various  captured  true  target  amplitudes  and  the  false  target  amplitude  as 
a  function  of  the  carrier  rf  phases  and  with  various  signal  magnitude  ratios  as  parame¬ 
ters.  R  is  seen  in  this  section  that  the  false  target  may  be  as  strong  as  the  true  targets. 
If  there  are  two  equally  strong  true  targets  and  one  smaller  true  target,  there  will  be  a 
false  target  of  approximately  the  same  size  as  the  smaller  true  target.  The  location  of 
the  false  target  changes  erratically  if  the  true  target  geometry  changes  slightly. 


APPROACH 

R  is  assumed  that  the  radar  transmits  a  phase-coded  signal  of  large  time-bandwidth 
product.  The  phase  codes  considered  In  this  report  are  linear  FM  and  linear  shift- 
register-generator  sequences,  as  they  are  described  for  example  in  Ref.  1.  The  target 
space  contains  a  number  of  discrete  point  scatterers  at  different  ranges;  that  is,  the  tar¬ 
gets  are  assumed  to  be  far  enough  separated  that  they  can  be  resolved  individually  by  the 
radar.  At  the  radar  receiver  input  there  will  therefore  be  a  summation  of  phase- 
modulated  sine  waves,  with  the  phases  between  the  sine  waves  depending  on  the  very  ac¬ 
curate  range  increments  between  the  multiple  targets  and  with  the  time  delays  between 
the  modulation  functions  depending  somewhat  less  sensitively  on  the  geometry.  The 
summation  of  the  individual  radar  returns  will  hence  be  both  amplitude  and  phase  modu¬ 
lated,  even  if  the  transmitted  signal  envelope  was  constant. 

The  type  of  transmitted  signal  calls  for  a  matched  fitter  or  a  correlator  as  a 
receiver-processor.  In  practice  the  receiver-processor  is  frequently  preceded  by  a 
hard-limiting  device,  which  may  be  operating  either  at  IF  or  on  the  in-phase  and  quadra¬ 
ture  components  of  bipolar  video  signals,  depending  on  the  radar  design.  Examples  of 
such  receiver  designs  are  given  in  the  next  paragraph.  For  the  purpose  of  this  analysis 
the  hard  limiter  is  assumed  to  be  at  IF.  The  hard  limiter  at  bipolar  video  can  be  handled 
as  a  special  case  of  the  IF  limiter,  wherein  the  input  signals  are  allowed  to  be  in  phase 
(positive)  or  180  degrees  out  of  phase  (negative)  and  are  not  allowed  to  have  phase  values 
in  between. 

To  keep  the  theoretical  model  as  simple  as  possible  without  losing  significance  it  is 
assumed  that  two  or  three  signals  with  various  relative  magnitudes,  delays,  and  RF 
phases  are  present  at  the  limiter  input,  ft  is  also  assumed  that  the  noise  is  negligible  at 
the  limiter  input.  The  assumption  of  a  large  signal-to-noise  ratio  may  not  always  be 
fulfilled,  and  in  such  cases  the  results  of  this  analysis  should  not  be  applied.  It  is  well 
known  that  the  limiter  acts  like  a  linear  device  causing  a  loss  of  only  1  to  1.5  dB  in  radar 
sensitivity  as  long  as  the  signals  are  sufficiently  deep  in  the  noise  (2,3).  This  study  is 
concerned  with  the  case  that  the  limiter  output  signal  can  be  predicted  from  the  radar 
and  the  target  parameters.  In  this  deterministic  case  the  nonlinearities  of  the  channel 
cause  the  formation  of  coherent  beat  products  which  may  correlate  with  the  transmitted 
radar  code  at  a  time  shift  which  does  not  correspond  to  the  actual  location  of  a  physical 
target.  In  this  case  a  spurious  target  response  is  generated.  The  amplitude  of  the 
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spurious  target  cannot  be  explained  as  simply  an  addition  of  range  side  lobes,  as  would 
be  indicated  if  Hnear-matched-filter  theory  would  hold.  The  false  targets  may  be  of  the 
same  magnitude  or  stronger  than  the  true  targets.  This  fairly  quantitative  claim  has 
been  supported  by  the  results  of  a  computer  simulation. 

Complex  signal  notation  is  used  throughout  this  analysis.  The  limiter  is  mathemati¬ 
cally  described  as  a  device  which  removes  the  amplitude  variation  from  the  complex 
signal.  The  limiter  output  is  a  complex  phase -modulated  signal  of  uniform  amplitude. 

It  is  assumed  that  the  limiter  operates  distort  ion -free,  i.e.,  that  the  phase  modulation  of 
the  input  signal  arrives  undistorted  at  the  output.  Much  of  the  approach  to  the  problem 
was  influenced  by  thoughts  presented  by  Nolen  in  a  paper  entitled  "Effects  of  Limiting  on 
Multiple  Signals"  (4).  To  tu  ^ .  t.  continuity  and  also  since  Nolen's  paper  is  not  gener¬ 
ally  available,  some  of  his  results,  particularly  those  pertaining  to  linear  FM,  are  re¬ 
viewed  in  this  report. 

It  may  be  argued  that  the  signal-to-noise  ratio  of  the  unprocessed  signal  is  ordinar¬ 
ily  very  small  in  typical  pulse-compression  or  phase-coded  CW  correlation  radar  sys¬ 
tems  of  large  correlation  gain.  Therefore  the  assumption  that  the  signal-to-noise  ratio 
is  large  would  in  many  cases  not  be  valid  and  the  conclusions  from  this  study  would  not 
apply.  It  is  true  that  the  signal-to-noise  ratio  may  be  exceedingly  small  for  a  minimum 
de  tectable  signal.  One  should  keep  in  mind,  however,  that  large  interfering  scatterers 
(clutter)  may  be  strong  enough  to  dominate  over  the  noise  at  the  limiter  input  even  if  the 
radsu  Jr  designed  to  detect  very  small  signals  by  virtue  of  a  high  correlation  gain. 


IMF LE MENTATION  EXAMPLES  OF  HARD-LIMITING 
RECEIVER-PROCESSORS 

Hard  limiting  may  take  place  at  IF  or  at  bipolar  video,  which  may  also  be  considered 
as  zero  IF.  Components  of  limited  dynamic  range,  like  RF  amplifiers  or  mixers,  may 
act  very  much  like  hard  limiters  as  soon  as  they  become  saturated  by  large  signals.  The 
IF  limiter  may  be  considered  as  a  device  which  ideally  would  preserve  the  phase  and  de¬ 
stroy  the  amplitude  modulation  of  a  signal.  If  the  input  signal  is  mathematically  de¬ 
scribed  by  a  complex  vector  of  variable  amplitude  and  variable  rotation  rate,  then  the 
output  signal  would  be  given  by  a  constant -amplitude  vector  which  points  at  any  given 
moment  in  the  same  direction  as  the  input  vector.  The  output  vector  may  be  assumed  to 
have  unit  amplitude  all  the  time.  To  practically  implement  a  hard-limiting  device  one 
may  use  amplifier  chains  whose  gain  is  controlled  through  a  feedback  loop  in  such  a 
fashion  that  it  is  inversely  proportional  to  the  input  signal  amplitude.  At  the  output  of 
the  amplifier  chain  one  would  observe  a  phase-modulated  sinusoidal  signal  with  constant 
amplitude  as  long  as  the  feedback  loop  is  fast  enough  to  respond  to  changes  of  the  input 
amplitude.  A  probably  less  troublesome  way  to  achieve  hard  limiting  is  to  use  the  cutoff 
characteristics  of  suitable  nonlinear  elements  like  transistors  or  diodes,  after  the  input 
signal  has  been  sufficiently  preampUfled. 

Figure  1  is  a  simplified  diagram  of  a  linear  FM  pulse-compression  receiver.  The 
RF  signal  coming  from  the  duplexer  Is  preampUfled,  heterodyned  to  a  conveniently  se¬ 
lected  IF,  pasted  through  s  hard  limiter  and  a  weighting  filter,  pulse  compressed  in  an 
ultrasonic  dispersive  delsy  line,  detected,  and  video  amplified  for  display.  The  weight- 
inr  filter  may  serve  a  dual  purpose  in  this  case.  B  may  he  used  to  reduce  the  range  aide 
lobes  of  the  compressed  radar  signal  and  to  eliminate  higher  harmonics  which  may  be 
generated  by  the  diode  limiter.  In  other  words  it  smooths  off  the  corners  of  the  signal 
coming  from  the  ttmiter.  A  typical  reason  for  placing  a  Hmttor  ahead  of  the  pulse- 
compression  device  la  to  normalise  the  noise  power  to  that  It  is  possible  to  set,  after 
pulse  compression  sod  rectification,  a  detection  threshold  at  a  computed  level  to  achieve 
a  constant  false  alarm  rate,  or  CFAR,  no  matter  how  strong  the  noise  or  interference  is 
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Fig.  1  -  Typical  linear  FM  pulse-compression  receiver 


before  limiting.  The  dynamic  range  of  the  radar  signal  becomes  increased  through  pulse 
compression  by  as  much  as  its  correlation  gain.  Limiting  may  also  be  used  to  keep  the 
dynamic  range  of  the  processed  signal  between  convenient  boundaries.  Another  reason 
for  limiting  may  be  that  it  is  simpler  or  cheaper  to  use  components  with  a  relatively 
small  dynamic  range  and  that  any  large  dynamic  range  after  pulse  compression  would 
not  be  needed  anyway. 

Figure  2  shows  a  commonly  implemented  layout  for  a  digital  correlator  using  binary 
shift  registers  as  memory  elements  and  operating  on  coherently  detected,  so  called  bi¬ 
polar  video  signals,  the  in-phase  or  I-signal  and  the  quadrature  or  Q -signal.  To  reduce 
equipment  complexity  one  may  omit  the  Q -signal  channel.  One  loses  on  that  case,  how¬ 
ever,  on  the  average,  3  dB  in  radar  sensitivity.  This  type  of  correlator  is  frequently 
used  to  process  pseudo -randomly  zero-pi -phase-coded  radar  signals.  One  of  its  ad¬ 
vantages  is  its  flexibility,  since  the  code  memory  may  contain  virtually  any  sequence  of 
plus  a:xl  minus  bits.  The  only  information  that  is  recorded  in  the  shift  registers  about 
the  rndar  signals  is  the  polarity  of  the  bipolar  video  signals  at  the  instant  of  sampling. 
The  value  of  its  amplitude  is  disregarded.  The  output  signals  are  therefore  the  same  as 
it  hard  limiting  had  taken  place  In  the  video  amplifiers  between  the  coherent  detector 
and  the  shift  registers  in  Fig.  2. 


GENERAL  THEORY  OF  THE  HARD  LIMITING  OF  TWO 
INPUT  SIGNALS 

The  approach  in  the  case  of  hard  limiting  of  two  Input  signals  is  the  same  as  the  one 
selected  by  Nolen  (4).  In  addition  to  reviewing  the  cases  of  two  constant -frequency 
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sinusoids  and  two  linearly  swept  FM  signals  which  have  been  treated  by  Nolen,  a  method 
of  handling  the  case  of  the  pseudo-randomly  zero-pi-phase-coded  signals  by  essentially 
the  same  technique  will  be  shown. 

Two  signals  at  the  limiter  input  may  be  described  mathematically  as  the  summation 
of  two  complex  vectors.  Through  separating  out  a  factor  .  xP  <  /  jtt /,,  r , ,  where  f„  may 
be  called  the  carrier  frequency,  one  displays  only  the  variations  of  the  complex  vectors 
with  respect  to  an  average  position  or  with  respect  to  a  reference  vector.  In  Fig.  3  there 
is  shown  as  an  example  two  superimposed  signal?  oi  different  magnitudes. 


Fig.  3  •  Phaaor  diagism  illustrating 
limiter  operation  in  the  complex  plane 


The  signal  with  the  burger  amplitude  will  be  arbitrarily  designated  signal  1.  If  the 
vectors  represent  two  sine  waves  of  constant  but  different  frequencies,  with  /  being  the 
frequency  of  sine  wave  1,  vector  1  may  be  considered  to  be  fixed  and  vector  2  to  be  ro¬ 
tating  with  Ute  difference  frequency.  The  true  (real)  electrical  signal  may  be  visualised 
ss  the  projection  on  the  real  axis  of  the  complex  vector  summation  rotated  around  the 
origin  at  the  rate  i„. 

The  heavily  drawn  vector  in  Fig.  3  represents  symbolically  the  Umiter  output  signal. 
Its  endpoint  always  falls  on  the  unit  circle  around  the  origin,  sod  it  is  aligned  in  parallel 
with  the  complex  vector  resultant  from  the  Unear  combination  of  the  input  signals. 

To  solve  the  problem  one  has  to  represent  the  Umiter  output  signal  aa  a  function  of 
the  phase  difference  of  Ute  input  signals  with  the  ratio  of  the  small  signal  amplitude  to 
the  large  signal  ampUtude  ss  s  parameter.  The  evolving  nonlinear  expression  is  too 
complicated  to  be  directly  useful,  however.  Developing  the  relation  between  the  output 
signal  and  the  phase  change  of  the  input  signals  in  a  Fourier  series  permits  a  much  more 
useful  functional  presentation  of  the  output  signal  In  the  foi.n  of  superimposed  coherent 
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phase-modulated  signals  with  amplitudes  that  can  be  calculated  essentially  through  eval¬ 
uation  of  Fourier  coefficients. 

Let  it  be  assumed  that  the  limiter  input  signals  have  the  constant  amplitudes  <l0  and 
4,  and  the  variable  phases  an  ’  o,(f),  vhere  t  is  the  time.  The  input  signal  sln 

which  is  the  sum  of  the  two  sign*' is  would  exhibit  both  amplitude  and  phase  changes  as  a 
function  of  time.  By  convention  the  amplitude  4,  shall  be  no  larger  than  40.  One  can 
then  define  the  small  signal  to  large  signal  amplitude  ratio  a,  which  would  never  be 
larger  than  1: 

,4 

a  r.  -  .L  <  1  (1) 


The  input  signal  may  then  be  expressed  as 


e  1  f  i4.  e  1 


’in  --  °  *  Al 


."V  *>  R  cja(t) 


In  the  last  line  of  this  formula,  R  represents  the  length  of  the  resultant  vector  described 
by  the  terms  between  braces  on  the  second  line,  and  a  represents  its  phase.  Reference 
is  made  to  Fig.  4  to  explain  the  relationship.  To  the  end  of  a  unit  vector  parallel  to  the 
real  axis  is  attached  a  smaller  vector  of  length  a  and  at  an  angle 


fl(f)  <!,,(' 


Fig.  4  -  Phasor  relationships 
in  the  complex  plane 
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The  action  of  the  ideal  limiting  device  is  to  replace  the  time-variable  amplitude  of 
the  signal  by  a  unit  amplitude  and  to  leave  its  phase  untouched.  This  may  be  mathemati¬ 
cally  accomplished  by  a  simple  omission  of  the  factors  40  and  R  in  Eq.  (2).  The  limiter 
output  hence  is  given  by 


>i*0<  *>♦“{*>] 

sout  ~  c 


(4) 


The  function  a  depends  on  the  phase  difference  e,  which  in  turn  is  a  function  of  t . 
From  the  geometry  depicted  in  Fig.  4  or  through  evaluation  of  the  identity 


oi  v  ;  obtains 


K  =  1  +  a  e'9  , 


.  a  sin  0 

a  -  tan  1  - - 

1  +  a  cos  0 


(5) 

(6) 


A  plot  of  a  as  a  function  of  o  and  with  the  parameter  a  ranging  from  0  to  1  in  steps  of 
0.2  is  pictured  in  Fig.  5.  If  a  -  o,  then  a  is  identically  equal  to  0.  For  small  values  of 
a,  say  for  a  =  0.2  the  function  resembles  a  sine  wave.  For  a  =  l,  a  is  a  linear  saw¬ 
tooth  function  connecting  from  a  =  -90  degrees  to  +90  degrees  and  with  the  discontinuity 
at  e  -  +180  degrees.  For  values  of  a  in  die  range  between  0.2  and  1  the  function  a (4) 
resembles  a  distorted  sine  wave. 


ALPHA  (DEGREES) 


Fig.  5  -  Relation  between  the  phase  of  the  sum  signal  and  the 
phase  difference  of  the  component  signals 


The  phase  difference  o  may  change  according  to  a  pseudo-random  sequence  as  a 
function  of  t,  or  it  may  be  some  other  very  complicated  sequence.  Hence  the  nonlinear 
expression  which  one  obtains  for  the  output  signal  through  combining  Eqs.  (4)  and  (6), 
namely, 


S 


out 


>[ 


exP  j  U0(f> +  t80'1 


a  sin  0(f) 

1  +  a  cos  0(f), 


(7) 


does  not  directly  indicate  which  signal  components  are  present  in  sout.  In  particular  it 
does  not  show  how  strongly  the  original  signals  are  present  and  whether  and  to  what 
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extent  new  signals  are  generated.  It  Is  a  fortunate  circumstance  that  the  Fourier  series 
development  of  Eq.  (7)  leads  to  a  summation  whose  terms  may  be  identified  in  several 
cases  of  practical  importance  with  images  of  the  input  signals  and  with  newly  generated 
signals.  Since  a  is  a  periodic  function  of  e  with  the  period  2v,  one  may  use  the  devel¬ 
opment 


=  exp  (j  tan"1 


a  sin  0  ' 

1  +  a  cos  9  , 


=  £  c„<a) 


-CD 


eJn 9 


(8) 


This  formula  expresses  an  identity  except  at  those  points  where  the  function  on  the  left 
side  has  a  discontinuity.  A  finite  number  of  summation  terms  may  provide  a  very  good 
approximation  except  in  the  vicinity  of  discontinuities.  It  may  be  remarked  that  Eq.  (8) 
does  not  represent  a  spectral  decomposition  of  the  limiter  output  but  rather  a  series  ex¬ 
pansion  for  a  nonlinear  relationship. 

All  coefficients  <^(a)  are  real,  since  a  is  an  odd  function  of  9.  The  proof  is  as 
follows:  As  a  consequence  of  the  relationship  ,  where  the  asterisk 

indicates  complex  conjugate,  one  may  equate  2  C„(a)  e~'n*  and  [£C„(a)  ein^]* .  Hence 
cn(a)  =  cji(a) ;  i.e.,  the  coefficients  C„  are  real.  In  the  general  case  the  coefficients  c„ 
and  c_n  will  not  be  the  same,  however. 

Nolen  has  shown  how  the  coefficients  C„  may  be  obtained  by  collection  of  terms 
from  an  infinite  product  of  infinite  series  in  powers  of  exp  (j9).  It  is  possible  to  calcu¬ 
late  the  C„  values  through  numerical  integration  methods,  which  may  perhaps  be  more 
easily  adapted  to  automatic  computer  evaluation.  The  coefficients  may  be  obtained  in  the 
usual  way  through  multiplication  of  Eq.  (8)  with  exp  {-jm9)  and  integration  over  9 : 


Hence 


♦  CP 

T\  C„(«)  e*1*  e-"*  d9 

fl*-CP 


ej(n-m)9  d6 


=  2vCm  ( a)  . 


c„(«) 


e/[a<0)-n9]  de  , 


This  integral  may  be  decomposed  into  a  real  and  an  imaginary  part: 


(9) 


(10) 


cos  [a(0)  -  nO]  de  +  j 


&  sin  [8(9) -nS] 

J -tr 


(11) 
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The  second  integral  is  equal  to  0,  since  a  is  an  odd  function  of  e.  By  the  same  token 
the  remaining  integral  taken  from  o  =  -v  to  o  =  +w  is  equal  to  twice  the  integral  taken 
from  e  =  o  to  o  -  v\ 


:  =  J-T 

"  2  TtJ., 


cos  [a (0)  -  n 6]  ad 


=  |  cos  [a(0)  -  no]  do  .  (12) 

Jo 

Inserting  Eq.  (6)  for  a  as  a  function  of  o  leads  to  the  full  expression  for  the  c„  values: 

c.,lf  cos  (tan-*  ■■  a  i»n-g-  -  n$)  de  .  (13) 

n  71  J0  \  1  ♦  a  cos  0  I 


a  cos  0 

Although  this  integral  cannot  be  evaluated  in  closed  form  for  any  arbitrary  value  of  the  sig¬ 
nal  amplitude  ratio  a,  it  can  be  integrated  in  the  special  case  a  =  i.  In  this  case  one 
obtains 

2  sin  4  cos  4 
sin  6  .  .  i  2  2 

a  (0)  =  tan  1  - =  tan  1 - 


=  tan-1 


1  +  cos  0 


6 

,,nl_  0 


2  cos*  4 


cos  ~2 


0  "  2 


(14) 


This  is  the  equation  for  the  straight  line  which  appears  in  Fig.  5  if  the  parameter  a  is 
equal  to  1.  Inserting  Eq.  (14)  into  Eq.  (12)  leads  to 


Cn(1)  =  *  1  cos(i  ■ n9) de J-f 


-Dn 


2n 


(15) 


Also  if  a  .  o,  i.e.,  if  the  smaller  signal  disappears  entirely,  a  trivial  solution  exists.  In 
that  case  a  is  identically  equal  to  0  and  one  obtains 


rn 

Cn(°>  =  | 

**0 


cos  no  do  -  1  .  if  n  =  0  , 


=  0,  for  all  other  n  values 


(18) 


The  general  shape  of  the  curves  in  Fig.  5  suggests  that  a  sine  function  might  be  a  reason¬ 
ably  good  approximation  for  a(0)  as  long  as  a  is  small  enough.  Analytically  one  may 
derive  from  Eq.  (6)  that  a  good  approximation  is 

a  =  a  sin  o  .  (17) 

for  small  values  of  a.  Inserting  this  approximation  into  the  integral  for  C„  leads  to 


if 

J  0 


cos  (a  sin  0  ~  nB)  d$ 


(18) 
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This  is  an  integral  representation  of  a  Bessel  function  of  the  first  kind  of  the  order  n  and 
with  the  argument  *  (Ref.  5). 

A  good  approximation  to  Ct  and  c_,  for  sufficiently  small  values  of  a  is  therefore 
C,  =  y,  (a)  and  c_,  =  yM(a).  From  the  first  term  of  the  McLaurin  series  for  J (x)  one 
obtains  the  approximate  relationships 

=  -  a/2  and  =  a/2  ,  if  a  <<  1  .  (19) 

It  would  be  erroneous,  however,  to  regard  Eq.  (18)  as  an  approximation  for  the  higher 
order  coefficients  C„.  The  higher  order  coefficients  C„  depend  very  critically  on  the 
higher  order  terms  in  the  development  of  a  as  a  function  of  e ,  and  exactly  those  have 
been  neglected  in  the  approximation  given  as  Eq.  (17).  A  valid  approximation  may  be  ex¬ 
pressed  as  a  summation  of  products  of  Bessel  functions  (4). 

The  approximate  expression  given  as  Eq.  (19)  along  with  the  special  result  given  as 
Eq.  (15)  provides  some  insight  into  the  general  behavior  of  the  Fourier  coefficients  c„ 
as  functions  of  a.  A  numerical  method  based  on  a  fast  Fourier  transform  computer 
program  has  been  used  to  obtain  numerical  answers  for  a  set  of  different  parameters. 

The  computer  program  and  the  methods  by  which  it  was  checked  are  explained  in  the  Ap¬ 
pendix  of  this  report.  The  results  are  presented  as  a  set  of  curves  in  Fig.  6  and  in  the 
form  of  line  spectra  in  Fig.  7.  Figure  7  also  exhibits  digital  printouts  for  the  c„  values. 
The  following  section  of  this  report  will  illuminate  the  physical  significance  of  the  Fourier 
coefficients  cn  in  a  special  case. 


Fig,  6  -  The  Fourier  coefficients  c„  as  functions  of  the  sig¬ 
nal  intensity  ratio  •  (labeled  A  by  the  computer)  and  of  the 
order  n  (labeled  N)  ranging  between  -4  and  5 
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Fig.  7(b)  -  Spectral  line  presentation  of  the  Fourier  coefficients  c  with  the  order  n  between  -14  and  15  and  strength  above  -50  dB 


-2  DB 


Fig.  7(d)  -  Spectral  line  presentation  of  the  Fourier  coefficients  cn  with  the  order 
n  between  -14  and  15  and  strength  above  -50  dB 
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Pig.  7(e)  -  Spectral  line  presentation  of  the  Fourier  coefficients  c„ 
with  the  order  n  between  -14  and  15  and  strength  above  -50  dB 
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Fig-  7(1)  -  Spectral  line  presentation  of  the  Fourier  coefficients  c. 
with  the  order  »  between  -14  and  15  and  strength  above  -50  dB 


An  expression  fer  the  limiter  output  signal  which  will  later  prove  to  be  quite  useful 
will  be  derived.  The  combination  of  Eqs.  (3),  (4),  and  (8)  results  in 


.♦to 


(20) 


where  the  coefficients  c„  are  the  Fourier  coefficients  which  have  been  discussed  in  this 
section. 


LIMITING  OF  A  PAIR  OF  CONSTANT-FREQUENCY  CW  SINUSOIDS 

The  phases  of  two  input  signals  that  are  constant-frequency  CW  sinusoids  can  be 
expressed  through  two  linear  functions  of  tins: 


*>,(»)  2 *y  «nri  *,(l)  2» /,#  .  *l#  .  (21) 

where  f.  and  are  the  constant  frequencies  of  the  sinusoids,  is  the  phase  shift 
which  exists  at  the  moment  t  o,  sad  the  small  to  larger  signal  amplitude  ratio  is 
given  by  the  parameter  *  which  has  been  defined  in  Eq.  (1). 

Insertion  of  Eq.  (21)  into  Eq.  (20)  permits  one  to  write  sa  expression  tor  the  limiter 
output: 


^♦(*"1  %»•*( r,  * **»••*»• » 


(22) 


This  is  aa  expression  tor  a  linear  superposition  of  constant -frequency  sinusoids,  consist¬ 
ing  of  frequency  /,  with  amplitude  c,<*),  frequency  ft  with  amplitude  c,<*),  and  side- 
bends  arranged  on  both  sidee  of  the  pair  of  input  frequencies  at  regular  intervals  deter¬ 
mined  by  the  frequency  difference  fx  r,  and  with  amplitudes  tapering  off  and  given  by 
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on  the  side  adjacent  to  /,  and  given  by  c_,, 


on  the  other 


side.  Consequently  Fig.  7  may  be  regarded  as  line  spectra  of  the  limiter  output  wave¬ 
forms  in  the  case  of  two  monochromatic  input  oscillations  with  different  frequencies. 

The  figure  shows  how  the  spectral  envelope  becomes  distorted  unsymmetrically  unless 
the  two  input  signals  are  equally  intense.  It  should  also  be  noted  from  Fig.  7  that  the 
higher  order  beat  products  (n  farther  from  0)  fall  off  much  faster  for  small  .i  values 
(s,  *0  closer  to  0  db)  than  for  a  values  in  the  vicinity  of  1.  It  is  evident  that  the  smaller 
signal  may  be  suppressed  relative  to  the  larger  signal  by  as  much  as  approximately  6  dB, 
as  it  is  well  known  (6). 

It  can  be  seen  directly  from  Eq.  (19)  that  the  amplitude  ratio  of  the  smaller  and  the 
much  larger  signal  after  limiting  is  equal  to 


-  —  .  «itii  .i  i 


(23) 


That  is  1/2  the  amplitude  ratio  before  limiting.  The  factor  1/2  in  amplitude  corresponds 
to  -6.02  dB.  The  smaller  signal  is  accompanied  by  a  mirror  image  of  approximately  the 
same  size. 


1 

m. 
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It  is  also  possible  to  explain  the  6-dB  suppression  of  a  snail  signal  without  going 
into  mathematics.  The  small  signal  nay  be  visualized  as  a  snail  vector  attached  to  the 
endpoint  of  the  larger  signal  vector.  The  vector  resultant  from  the  summation  of  both 
input  vectors  equals  the  large  signal  vector  with  a  small  additional  amplitude  and  phase 
modulation  caused  by  the  smaller  signal.  Half  of  the  small  signal  energy  produces  the 
amplitude  modulation  and  half  produces  the  phase  modulation  of  the  signal  superposition. 
After  limiting,  the  signal  energy  is  normalized  and  the  amplitude  modulation  is  removed, 
which  would  account  for  a  3-dB  loss.  The  phase  modulation  is  contained  in  two  sidebands, 
namely  in  the  captured  small  signal  and  in  its  image.  Both  are  of  approximately  the 
same  intensity  and  thus  suffer  a  power  splitting,  which  accounts  for  another  3-dB  loss  in 
the  power  balance. 

A  simple  bench  test  has  been  run  by  R.  M.  Crisler  to  check  on  the  validity  of  the 
theoretical  results  described  in  this  section.  Two  constant-frequency  sine  waves  of  high 
spectral  purity,  one  at  60  MHz  and  the  other  one  at  60.016  MHz  were  added  together  in  a 
hybrid  circuit,  amplified  in  an  intentionally  strongly  overdriven  and  hence  hard-limiting 
IF  strip  and  displayed  and  photographed  on  a  Singer  Model  SPA-100  microwave  spectrum 
analyzer.  The  amplitude  ratio  «  of  the  two  sinusoids  could  be  adjusted  by  means  of  a 
calibrated  attenuator  inserted  in  the  60.016-MHz  signal  path  before  the  hybrid.  The  IF 
amplifier  may  be  assumed  to  have  a  fairly  flat  amplitude  response  over  the  frequency  in¬ 
terval  of  interest  Figure  8(a)  shows  the  signal  superposition  before  entering  the  IF 
strip  as  it  appeared  on  the  screen  of  a  Tektronix  type  454  oscilloscope.  In  Fig.  8(b)  is  a 
display  of  the  amplifier  output,  which  exhibits  fairly  constant  amplitude  except  in  the 
Immediate  vicinity  of  the  mills  of  the  input  signal  envelope.  When  the  input  signals  dif¬ 
fered  by  a  few  decibels,  the  minima  of  the  unlimited  signal  envelope  were  large  enough  to 
saturate  the  IF  amplifier  so  that  one  would  observe  a  constant  amplitude  oscillation  with¬ 
out  any  dips  In  the  output  Figures  8(c)  through  80)  show  the  results  of  the  spectral  anal¬ 
ysis  of  the  linearly  superimposed  sinusoids  (Fig.  8(c))  and  of  the  hard-limited  signals 
(the  rest  of  the  pictures).  The  vertical  scale  is  linear  In  decibels,  each  division  corte- 
spomttng  to  a  10-dB  step.  Figure  8(c)  shows  that  the  input  frequencies  are  free  from 
spurious  sidebands  before  they  enter  the  limiter,  fin  Fig.  8(d)  one  sees  s  symmetrical 
array  at  spectral  Haas  generated  oa  both  sides  of  the  equally  strong  input  frequencies, 
ft  should  be  remarked  that  the  amplitudes  of  the  signal  frequencies  t%  sad  f.  ere  ap¬ 
proximately  4  dft  below  the  intensity  provided  by  the  limiter  for  a  single  input  signal 
(compare  Fig.  8(d)  with  Fig.  8Q)). 
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(a)  Scope  picture  of  two  equally 
strong  sinusoidal  signals  be¬ 
fore  limiting 


(e)  Spectrum  of  the  limiter 
output  signal  from  two  input 
signals  differing  by  5  dB 
(same  as  Fig.  7(d)) 


(h)  Spectrum  of  the  limiter 
output  signal  from  two  input 
signals  differing  by  20  dB 
(same  as  Fig.  7(g)) 


(b)  Scope  picture  of  the  same 
signals  as  in  Fig.  8(a)  but  after 
hard  limiting 


(f)  Spectrum  ef  the  limiter 
output  signal  from  two  input 
signals  differing  by  10  dB 
(same  as  Fig.  7(e)) 


(I)  Spectrum  of  the  limiter 
output  signal  from  two  input 
signals  differing  by  28  dB 
(name  ae  Fig.  7(h)) 


(c)  Spectrum  of  the 
signals  of  Fig.  8fs) 


(g)  Spectrum  of  the  limiter 
output  signal  from  two  input 
signals  differing  by  IS  dB 
(Mine  as  Fig.  7(1)) 


drum  of  the  limiter 
signal  from  one  iaput 


(d)  spectrum  of  the  limiter 
output  rtgnal  (Fig.  8(b))  from 
two  equally  strong  input  Mg- 
nals  (same  as  Ftg.  7(a)) 


Fig.  S  -  Waveforms  sad 


recorded  is  a  beech  tot 
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If  the  input  signals  differ  by  as  little  as  5  dB,  the  stronger  output  signal  obtains 
nearly  its  full  size,  and  the  smaller  signal  appears  at  approximately  -10  dB.  That  means 
the  smaller  3ignal  was  suppressed  by  approximately  5  dB  through  the  limiter  and  it  is 
accompanied  by  a  slightly  smaller  image  line  and  by  a  regularly  spaced  array  of  spectral 
lines  which  taper  off  unsymmetrically  (Fig.  8(e)).  if  the  input  signals  differ  by  10  dB 
(Fig.  8(f)),  then  nearly  the  full  amount  of  the  theoretically  predicted  6-dB  small  signal 
suppression  may  be  observed.  Figures  8(g),  8(h),  and  8(i)  show  the  small  signal  sup¬ 
pressed  by  6  dB  and  accompanied  by  an  equally  strong  image  when  the  input  signals  dif¬ 
fered  by  15,  20,  and  25  dB  respectively  before  limiting.  Figure  8(j)  shows  a  single  spec¬ 
tral  line  representing  the  fundamental  frequency  of  a  hard-limited  single -frequency  input. 
The  spectral  lines  displayed  in  Figs.  8(d)  to  8(i)  are  the  beat  frequencies  of  the  two  input 
signals  due  to  the  nonlinear  processing;  the  spectral  lines  are  not  harmonic  frequencies. 
Harmonic  frequencies  which  are  present  if  a  60-MHz  signal  is  limited  would  be  at  multi¬ 
ples  of  60  MHz,  arJ  they  would  therefore  be  outside  the  spectral  range  displayed,  which 
covered  an  interval  of  approximately  200  kHz  centered  around  60  MHz. 

Point-by-point  comparison  of  the  measured  spectra  such  as  shown  in  Fig.  8  and  the 
calculated  c„  values  pictured  in  Fig.  7  discloses  their  exact  matching. 

This  section  of  the  report  may  be  concluded  with  the  remark  that  the  beat  frequencies 
generated  through  hard  limiting  would  appear  like  false  targets  in  the  case  of  a  CW  radar 
or  a  pulse  doppler  radar  of  insufficient  dynamic  range.  As  applied  to  a  pulse  doppler 
radar  one  may  say  that  two  clutter  components  of  slightly  different  doppler  frequencies 
but  both  falling  in  the  clutter  notch  of  the  processor  response  —  if  strong  enough  to  drive 
the  IF  amplifier  into  saturation  —  may  generate  beat  frequencies  which  would  be  within 
the  accepted  doppler  domain. 

It  may  be  noted  also  that  an  analogy  exists  with  an  antenna  problem.  False  angular 
responses  (ghost  targets)  ar  e  observed  in  a  two-target  environment  if  the  elements  of  an 
array  antenna  are  nonlinear  (7,8). 


LIMITING  OF  A  PAIR  OF  LINEARLY-FREQUENCY- 
MODULATED  SIGNALS 

Consider  a  pair  of  linearly-frequency- modulated  signals  in  the  idealized  case  where 
the  input  signals  are  entirely  overlapping  and  very  strong  as  compared  to  thermal  noise. 
The  analysis  is  very  similar  to  the  one  given  in  the  previous  section  for  the  case  of 
constant-frequency  sinusoids.  A  pair  of  linear  FM  signals  gives  rise  to  beat  frequencies 
within  the  bandwidth  of  the  frequency  excursion.  These  constitute  signals  which  are 
compressed  by  the  matched  filter  or  correlator  and  which  will  appear,  after  processing, 
as  an  array  of  false  targets  both  behind  and  ahead  of  the  pair  of  true  targets  (4).  If  the 
expanded  signals  overlap  less  than  completely,  there  will  be  a  proportionate  reduction  of 
the  false  target  magnitudes  (4).  Let  one  input  signal  be  given  by 

A0  exp  j  (2tj  f0t  +  bt J) 

and  the  other  by  a  scaled,  time-  and  phase-shifted  version  of  the  first  signal: 

4,  exp  /  [2trf0{t-  T)  +  b(t-T)2  +  <i,0]  . 

The  Instantaneous  frequency  of  the  first  signal,  which  is  by  definition  equal  to  the  deriva¬ 
tive  of  its  phase  divided  by  2v,  is  then  equal  to 
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and  the  instantaneous  frequency  of  the  second  signal  is  given  by 

'o  ‘  W(‘-T >  • 

where  r  is  the  mutual  time  shift  between  the  two  input  signals.  Inserting 

<V * )  2 irf0t  *  ht 2 

and 

T  2 ^f^(l-T)  i  b  (t  -  T) 2  ♦  «i10 

into  Eq.  (20)  results  in 

c  (2wr0(»-nr)4h(<-nr)J«n<i10-(n-l)nftrJ] 


Equation  (26)  says  that  the  limiter  output  may  be  decomposed  into  a  summation  of  linear 
FM  signals.  By  comparison  with  the  expressions  for  the  input  signals  one  derives  that 
the  limiter  output  consists  of  an  array  of  mutually-time-shifted  images  of  the  input  wave¬ 
forms,  with  delay  r  between  adjacent  components  and  with  amplitude  c„  for  the  nth 
component. 

At  the  output  of  the  matched  filter  or  correlator  which  would  follow  in  the  signal 
processing  scheme  after  the  hard  limiter  one  would  therefore  observe  an  array  of  com¬ 
pressed  pulses,  namely,  the  large  signal  response  with  amplitude  c0,  the  smaller  signal 
response  with  amplitude  c,,  an  image  response  with  amplitude  c_,,  and  smaller  pulses 
at  regular  intervals  and  with  amplitudes  c2,  c,,  . . .  on  the  side  of  the  smaller  true  sig¬ 
nal  and  with  amplitudes  c_2 ,  c.3,  ...  on  the  side  of  the  image.  If  the  compressed  and 
detected  signals  were  passed  through  a  logarithmic  amplifier,  one  would  be  able  to  ob¬ 
serve  at  its  output  essentially  the  same  waveforms  as  the  ones  depicted  in  Figs.  8(d) 
through  8(j).  All  that  has  been  said  in  the  previous  section  on  the  amplitude  ratios  and 
capture  effects  would  apply  also  in  the  case  of  linear  FM.  One  could  in  fact  envisage  the 
spectrum  analyzer  displays  8(d)  through  8(j)  also  as  scope  pictures  for  the  case  of  linear 
FM,  simply  by  considering  the  traces  as  functions  of  time  instead  of  frequency.  It  is 
then  obvious  that  one  observes  false  targets  in  addition  to  the  true  targets.  This  phe¬ 
nomenon  may  be  very  cumbersome  if  the  radar  should  be  used  for  more  functions  than 
simple  target  detection.  The  way  out  of  this  dilemma  would  of  course  be  to  use  a  re¬ 
ceiver  and  signal  processor  of  sufficiently  large  dynamic  range. 

Attention  should  also  be  given  to  the  phase  term  2 vfnnT  in  Eq.  (26),  which  causes 
the  doppler  frequency  of  false  targets  to  increase  with  order  n  if  there  is  a  mutual 
motion  between  the  true  targets.  This  means  that  the  doppler  of  false  targets  may  be 
within  the  acceptance  domain  of  a  pulse-compression  pulse  doppler  radar  even  if  the 
"true"  responses  due  to  clutter  would  be  rejected  by  means  of  doppler  filtering  of  the 
compressed  signals. 


j 

1 


(24) 

(25) 

(26) 
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LIMITING  OF  A  PAIR  OF  ZERO-PI-PHASE-MODULATED  SIGNALS 

Consider  now  the  class  of  radar  signals  which  may  be  generated  through  switching 
the  phase  of  a  sinusoid  between  0  and  180  degrees  depending  on  a  binary  code.  The  bi¬ 
nary  code  may  be  visualized  as  a  sequence  of  zeros  and  ones  or  as  a  sequence  of  plus 
and  minus  signs  or  as  a  video  type  waveform  alternating  between  two  discrete  voltage 
levels.  Since  the  result  of  this  section  will  be  generally  valid  for  any  binary  coded  se¬ 
quence,  no  further  restriction  will  be  made  at  this  point  as  to  specific  code  classes.  An 
expression  shall  be  derived  for  the  small  and  the  large  signal  size  after  hard  limiting  as 
a  function  of  the  RF  phase  <t>i0  between  the  input  sinusoids  and  with  the  input  amplitude 
ratio  a  as  a  parameter.  As  before  in  this  report  the  RF  signals  are  assumed  to  be  com¬ 
pletely  overlapping  and  strong  as  compared  to  the  noise.  Equation  (20)  will  permit  one 
to  decompose  the  limiter  output  into  a  meaningful  sum  of  component  signals  which  may 
be  identified  with  the  input  signals.  Let  the  code  be  represented  by  a  sequence  of  coeffi¬ 
cients  ck  whose  values  are  either  0  or  1.  The  phase  modulation  impressed  on  a  radar 
signal  may  then  be  expressed  as 


Mn 


rect 


(27) 


where  the  summation  runs  over  all  k  values,  t  is  the  time,  t  is  the  duration  of  a  bit, 
and  rect  x  stands  for  Woodward's  rectangular  function  notation: 


rect  x  =  1  ,  for  x  in  the  interval  (-0.5.  0.5)  , 
-  0  ,  outside  the  same  interval  . 


(28) 


K  the  carrier  frequency  is  given  by  fQ,  then  the  phase  of  the  stronger  input  signal  may 
be  written  as 

*0<O  =  2"V  +  *e(*>  •  (29) 

and  the  phase  of  the  smaller  signal  may  be  derived  from  the  phase  of  the  stronger  signal 
through  introduction  of  a  time  delay  T  and  the  mutual  RF  phase  <a10: 

4,(0  =  2*V  ♦  *  tf,0  •  (30) 


Inserting  these  expressions  into  Eq.  (20)  yields 

I  WO  =  £  c"  eitn^(,‘r)'<B*1)0c(°tn^'>1  „'*v  .  (31) 

| 

r 

\  It  is  now  important  to  observe  that  both  <AC(  f )  and  -  r>  can  assume  only  either  one 

f  of  the  discrete  values  0  or  ».  Hence 


?  c  c  =  e  for  all  even  values  of  n, 

=  \  for  all  odd  values  of  n  . 


(32) 


Equation  (31)  may  therefore  be  transformed  into  two  separate  summations,  one  of  them 
combining  all  even-index  terms  and  the  other  collecting  all  odd-index  terms: 
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wo 


+  ® 

L 


;  2>)<iiS0\  ><*,c ( t  )  j2nf0t 


C  e 
L2n  e 


+  0?  \ 

c2  J  e'<2n4l>^ioj  eW->  e>2w/0» 


(33) 


From  this  expression  one  can  read  that  the  limiter  output  may  be  written  as  the  summa¬ 
tion  of  two  signals,  each  of  them  being  of  the  same  functional  shape  as  one  of  the  input 
signals,  and  with  the  complex  amplitudes  described  by  the  sums  of  exponentially  weighted 
Fourier  coefficients  of  even  or  odd  order  respectively: 


(Ec,e'W'»)  or  (I)  C2n+1  e''(2n+I),il°)  .  (34) 

These  expressions  are  consequently  properly  called  the  amplitude  of  the  stronger  signal 
and  the  amplitude  of  the  smaller  signal  after  limiting  has  taken  place.  Since  Eq.  (33)  is 
a  complete  description  of  the  limiter  output  signal  as  long  as  the  assumptions  made  are 
valid,  there  is  no  indication  of  false  targets. 

It  is  possible  to  transform  expressions  (34)  back  to  their  original  domain  in  a  gen¬ 
eral  fashion  so  that  explicit  knowledge  of  Fourier  coefficients  c„(a)  will  not  be  required 
for  numerical  evaluation  of  capture  effects.  From  Eq.  (8)  one  may  derive 


e'°W)  =  £  C„(a)  e'"W)  =  £  (-1)"  Cn(a)  c>n0  .  (35) 

Hence  one  obtains  through  addition  of  Eqs.  (35)  and  (8) 

e>o(0)  4  e'«W>  =  2  £cjn(a)  eJin6  .  (36) 

Through  subtraction  one  obtains 

,  2  £  Cjn+J(a)  .  (37) 

These  are  general  relations  exactly  of  the  form  as  needed  to  express  the  limiter  output 
signal  amplitudes  (34)  in  terms  of  the  function  exp  ja(0).  The  large  signal  amplitude 
after  limiting  is 


i  2  rn'i 


'2n 


10  - 


/2 


and  the  small  signal  amplitude  is 


£ 


'2n*  J 


j(  2nU)<d,0 


i 

c 


«(</>, 0  ) 


io(<A 


10 


+i/2 


(38) 


(39) 


Equations  (38)  and  (39)  have  been  programmed  for  automatic  computer  evaluation.  The 
results  are  depicted  in  Fig.  9  both  for  the  large  and  for  the  small  signal  distributions. 
There  is  a  horizontal  line  at  the  -3-dB  level  labeled  with  the  parameter  0  dB  indicating 
that  the  power  is  equally  split  between  the  output  signals  if  the  input  signals  are  equally 
strong.  All  curves  below  or  above  the  -3-dB  line  refer  to  the  smaller  or  stronger  signal 
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U  45  90  135  180 


D1C  IN  0EGREE5 

Fig.  9  -  Capture  of  a  pair  of  zero-pi  phase-coded  signals  as  a 
function  of  the  carrier  rf  phase,  with  the  signal  intensity  ratio 
as  the  parameter 


respectively.  The  curves  indicate  that  the  capture  effects  are  minimized  if  the  RF  phase 
between  the  input  signals  is  around  90  or  270  degrees.  If  the  RF  phase  amounts  to  zero 
or  180  degrees,  the  smaller  signal  will  be  suppressed  entirely,  at  least  in  the  absence  of 
noise  or  any  other  type  of  third  signal  component.  Noise  in  addition  to  the  two  input  sig¬ 
nals  will  permit  the  smaller  signal  to  determine  the  limiter  output  signal  with  some  fi¬ 
nite  probability  even  in  the  case  of  an  RF  phase  difference  of  0  or  180  degrees.  If  there 
is  a  correlated  third  signal  at  the  limiter  input,  there  may  be  a  false  target  generation 
effect,  as  will  be  shown  in  the  next  section.  It  should  be  mentioned  that  Eqs.  (38)  and 
(39)  may  be  used  to  derive  a  simple  geometric  construction  method  in  a  complex  plane 
to  describe  the  capture  effects  of  two  zero-pi-phase-coded  signals.  It  is  worth  under¬ 
scoring  that  Eqs.  (38)  and  (39)  properly  describe  the  capture  effects  no  matter  what 
particular  binary  phase  code  is  used. 

GENERAL  THEORY  OF  THE  HARD  LIMITING  OF  THREE 
INPUT  SIGNALS 

The  previous  sections  have  discussed  and  provided  mathematical  descriptions  of 
how  the  limiting  of  a  pair  of  linear  FM  signals  gives  rise  not  only  to  capture  effects  but 
also  to  spurious  target  generation,  whereas  any  kind  of  zero-pl-phase-coded  signal  pair  is 
subject  only  to  capture  and  not  to  false  target  generation.  The  absence  of  false  targets  in 
the  case  of  zero-pi  phase  codes,  although  in  agreement  with  general  knowledge  on  this 
subject,  may  seem  somewhat  baffling  from  a  theoretical  point  of  view,  in  particular  since 
certain  pseudo-random  sequences,  namely  linear  maximal  sequences,  which  are  fre¬ 
quently  used  as  radar  modulation  functions  and  which  are  easily  generated  by  means  of 
feedback  shift  registers,  exhibit  a  very  strong  structure  and  have  the  property  to  produce 
a  sequence  of  the  same  code  but  with  a  different  delay  if  combined  linearly.  This  prop¬ 
erty  justifies  the  adjective  "linear"  in  "linear  maximal  sequences,"  and  when  one  speaks 
here  of  a  linear  code  combination,  one  means  the  new  code  generated  if  two  codes  with  a 
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relative  delay  between  them  are  combined  according  to  the  ruin  of  Boolean  algebra 
on  a  bit-to-bit  basis.  This  kind  of  linear  combination  would  also  take  place  if  two  phase 
modulated  RF  signals  are  multiplied  or  beat  together,  since  the  multiplication  of  com¬ 
plex  signals  results  in  the  addition  of  the  exponents.  The  effect  of  signal  multiplication 
takes  place  at  any  nonlinear  element  in  the  signal  channel  and  especially  in  a  limiter.  It 
was  therefore  speculated  that  false  targets  might  be  generated,  if  not  in  the  case  of  two 
input  signals,  in  the  case  of  three  or  more  input  signals. 

The  theory  developed  is  a  simple  extension  of  the  analysis  made  for  the  two-input- 
signal  case  presented  earlier  in  this  report.  The  same  assumptions  are  made  through¬ 
out,  namely,  complete  expanded  signal  overlap  and  high  signal -to-noise  ratio.  The  input 
signal  to  the  limiter  may  then  be  written  as 

S.n  =  ^  e'^  +  A,  eJ*'(t)  +  A2 


=  e 


*»<*>(■ 


H  f 


,  _  _i<t> 


'<t>2  \ 


=  R  e'a(f)  . 


(40) 


The  convention  is  that  the  first  signal,  A0  e^o(  *>,  is  the  strongest  signal,  so  that  it  is 
meaningful  to  describe  the  input  signal  s,„  as  the  strongest  signal  modified  in  amplitude 
and  in  phase  through  the  presence  of  two  smaller  signals,  mathematically  expressed  in 
Eq.  (40)  through  R  and  a.  Set 


.  ^2 
ai  rT*  and  aa  ZT~ 

Ao 


(41) 


a,  and  -t3  being  the  two  ratios  of  a  smaller  signal  to  the  large  signal,  and  set 

°l  z  \  and  *  <*0  ' 


(42) 


The  ideal  limiter  normalizes  the  amplitude  of  the  input  signal,  Eq.  (40),  while  preserving 
its  phase  function  <£0(  t )  +  a(  t ) : 


“out 


/[<40( ')♦«(')] 


(43) 


where  a  is  the  phase  angle  of  the  complex  factor  between  large  parentheses  on  the  sec¬ 
ond  line  of  Eq.  (40): 


a(t)  ~  ton"1 


n,  sin  Sj  ♦  sin  0 a 
1  ^  flj  cos  0 j  *  cos  0 j 


(44) 


To  obtain  a  physically  meaningful  interpretation  of  the  limiter  output  signal  it  has  to  be 
decomposed  into  a  summation  of  signals.  In  the  previously  treated  case  of  a  pair  of  in¬ 
put  signals  it  proved  successful  to  develop  a  as  a  function  of  »  into  a  Fourier  series. 
Applying  the  same  technique  in  the  case  of  a  triplet  of  input  signals  leads  us  to  consider 
the  development 
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m  =  -cp 


;(m0  4  nO-  ) 
e  1  1 


(45) 


which  is  a  double  Fourier  series  for  a  as  a  function  of  o  j  and  o2.  Combining  Eq.  (43) 
with  Eqs.  (45)  and  (42)  leads  to  the  series  development  for  the  limiter  output: 


out 


LL 


i  +  (  l-m-n)cA0] 


(46) 


The  Fourier  coefficients  cmn  may  be  expressed  as  a  double  integral,  and  they  may  be 
determined  by  numerical  integration  methods.  In  the  case  of  practical  importance  which 
will  be  studied  in  the  next  section  one  obtains  general  expressions  for  limiter  output  sig¬ 
nal  amplitudes  which  may  be  transformed  back  to  the  original  domain,  so  that  any  ex¬ 
plicit  knowledge  of  Fourier  coefficients  will  not  be  required.  The  Fourier  series  ap¬ 
proach  permits  one  to  derive  relationships  in  terms  of  the  original  quantities. 


LIMITING  OF  A  TRIPLET  OF  LINEAR  CODED  SEQUENCES 

In  the  case  of  a  triplet  of  linear  coded  sequences  let  the  phase  modulation  function  of 
the  strongest  signal  be  given  by 


<MO  =  rePLTV  £  Ck  feCt  1  t*"  '  (47) 

*=  1 

where  t  is  the  time,  r  is  the  duration  of  a  bit,  cfc  =  0  or  1  depending  on  a  maximal  lin¬ 
ear  sequence,  l  is  the  code  length,  and  reptT  indicates  the  periodic  repetition  of  the 
code,  the  repetition  period  being  equal  to  Lr.  A  modulation  function  of  this  kind  is  usually 
generated  by  means  of  shift  registers  with  suitable  feedback  connections.  A  shift  register 
generator  consisting  of  s  stages  may  be  used  to  generate  a  code  of  length  L  =  2s  -  1. 

This  kind  of  code  generation  is  frequently  used  in  radar,  both  for  theoretical  and  practi¬ 
cal  reasons.  This  method  requires  little  hardware,  and  it  yields  a  code  with  a  nearly 
ideal  autocorrelation  function,  the  range  side  lobes  being  at  a  uniform  level  of  -20  log  L 
dB  below  the  main  peak.  The  low  and  uniform-range  side-lobe  level  is  a  consequence  of 
the  strong  structure  inherent  in  this  kind  of  pseudo-random  sequence,  namely,  of  the 
property  that  the  codes  added  bit  by  bit  to  a  delayed  version  of  the  same  code  under  the 
rules  of  Boolean  algebra  will  result  in  the  original  code  sequence  shifted  by  a  pseudo¬ 
random  number  of  bits.  It  will  be  seen  that  the  same  property  is  responsible  for  the  gen¬ 
eration  of  a  false  target  of  sometimes  larger  magnitude  than  any  one  of  the  suppressed 
true  target  returns. 

As  an  example  consider  the  shift  register  generator  depicted  in  Fig.  10.  At  its  out¬ 
put  one  would  observe  a  periodic  binary  code  sequence,  the  repetition  period  being  equal 
to  31  clock  pulse  intervals.  Since  the  shift  register  contains  five  stages,  the  feedback  is 
connected  as  shown  in  Fig.  10  to  ensure  maximum  code  length  and  L  =  2S-1  =  25-1  = 
31.  If  one  would  start  to  count  a  new  period  for  example  when  the  shift  register  is  loaded 
with  the  sequence  1,  0,  0,  0,  0,  then  as  a  consequence  of  the  operation  depicted  in  the  fig¬ 
ure  one  would  be  able  to  record  the  output 

1000010101110110001111100110100. 

If  this  code  word  is  added  to  the  same  word  but  shifted  to  the  right  by  one  bit,  namely, 

0100001010111011000111110011010, 
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FEEDBACK 
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Fig.  10  -  Maximum-code- length  shift  register  code  generator 


one  obtains  as  a  result  of  the  binary  addition  the  new  code  word 

1100011111001101001000010101110. 

The  remarkable  feature  of  linear  coded  sequences  is  that  the  code  word  obtained  is  not 
any  arbitrary  sequence  but  is  the  same  as  the  original  sequence  shifted  by  a  number  of 
bits,  in  this  case  shifted  by  18  bits  as  one  can  easily  verify  by  inspection  of  the  sequences. 
Binary  code  addition  takes  place  if  a  signal  of  the  type  exp  )fe(  t  >  is  multiplied  with  an¬ 
other  signal  exp  t-T),  T  being  the  relative  delay  between  the  two  signals,  since  the 
product  of  exp  Hdt)  and  exp  /*c(f-T)  is  exp>[«c(f)  +  *c(t-r)].  In  the  special  case 
of  linear  coded  sequences  one  has  then 

_  eA&e<»-r*) 


and  also 


eM.<  * )  V  e  *-r  >  (48) 

r  being  a  delay  which  depends  in  a  pseudo-random  fashion  on  the  shift  between  the  codes 
added  together. 

Now  apply  the  results  from  the  previous  section.  The  limiter  input  signal  is  mathe¬ 
matically  described  by  Eq.  (40),  where  the  exponent  of  the  strongest  signal  is 

>40(f> s  ♦  2»V)  ,  (49) 

the  exponent  of  the  second  signal  is 

/*,<*> 1  i  (*c(»  -  r,>  ♦  2 »V  ♦  4,01  .  (50) 

and  the  exponent  of  the  third  signal  is 

>*,<»)  *  ♦  2 ♦  *„)  .  (51) 

Here,  /,  is  the  RF  frequency,  <*10  and  <*]t  are  the  RF  phases  between  the  second  and 
the  first  and  Hie  third  and  Hie  first  signal  respectively,  and  r,  and  r,  are  the  mutual 
delays,  inserting  Eqs.  (49),  (50),  and  (51)  into  Eq.  (40)  results  in 
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•  ££c„„< 


^  e-/(m«n-|)^<.(»)«/m<4c(f-r|  )«>n^c(»-T2)  i2w/0(») 


(52) 


This  expression  may  be  significantly  simplified  if  one  realizes  that  exp  jk6c  is  equal  to 
exp  y<6c  or  to  1,  depending  on  whether  k  is  odd  or  even.  Using  also  the  property  given 
by  Eq.  (48)  one  can  write  the  following  truth  table  for  the  exponential  expression 

exp  [-;  (m  +  n  -  l)  <AC(  t  )  +  t  -  7,  )  +  jn<t,c(  t  -  Tj  )] : 


even  m 

odd  m 

even  n 

eJ6c(t) 

odd  n 

!«&,(<  )♦«*,(  »-r,  >*«5e.<  *-**2  »1 

One  can  now  write  Eq.  (52)  in  the  form  of  a  sum  of  four  terms  in  such  a  way  that  each 
index  of  a  Fourier  coefficient  cmn  is  even  or  odd  throughout  each  summation: 


6Ut 


H  C»*l.  2»(  T 

I 


cj*.  a*«i<*i  *a>  e 

* 


j[(J**!K4|0‘ai<Jj0! 

#(a*ol0*  (ai«ihiI0! 


>Ae<  »-7.  ) 

r  c  1  f  0 


t 


#«M***‘>  #a»V 

r  t 


(53) 


This  expression  shows  that  the  limiter  output  consists  of  four  signals,  three  of  them  be* 
ing  Identical  with  the  three  input  signals  except  for  an  amplitude  change  and  the  fourth 
signal  being  an  image  of  the  coded  sequence  at  some  pseudo-random  starting  position. 
The  amplitudes  of  all  four  signals  can  be  calculated  if  the  ratios  of  the  small  to  large 
signals  before  limiting  and  the  RF  phases  are  known,  without  requiring  specification  of 
the  particular  linear  code  used.  The  amplitudes  are  given  by  the  expressions  in  the 
large  enclosures  in  Eq.  (S3).  These  are  also  the  amplitudes  which  the  four  correlation 
peaks  of  the  signal  *ou,  will  have  after  crosscorrelation  processing  with  respect  to  the 
transmitted  radar  code  r ; .  ft  is  hence  possible  to  derive  quantitative  information 
about  capture  and  false  target  generation  effects  by  evaluating  the  sums  of  Fourier  co¬ 
efficients  as  Indicated  by  Eq.  (S3). 


As  has  been  indicated,  the  Fourier  series  approach  reveals  Its  usefulness  sad  power 
through  the  fact  that  the  summations  in  Eq.  (S3)  can  be  transformed  back  to  the  original 
domain,  thereby  enabling  numerical  discussion  of  results  without  a  **d  to  know  the  Fou¬ 
rier  coefficients.  Through  inserting  the  arguments  <*,  •  *  and/or  «2  •  in  place  of  a , 

and/or  «,  in  Eq.  (48),  one  derives  the  following  three  relations: 
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(M) 


,‘-W>  -  Z  X>1>" 


(55) 


and 


(5«) 


These  expressions  contain  the  same  Fourier  components  as  Eq.  (45),  however,  with  al¬ 
ternating  signs,  so  that  coefficients  with  even  or  odd  *  and  n  values  may  be  made  to  re¬ 
inforce  or  cancel  by  suitably  combining  Eqs.  (45),  (54),  (55),  and  (56).  Thereby  a  back 
transform  for  the  amplitude  relations  in  Eq.  (55)  is  obtained. 

In  adding  Eq.  (45)  and  Eq.  (54)  all  the  terms  with  even  *  values  will  be  doubled  and 
all  the  terms  with  odd  •  values  will  be  canceled.  This  may  be  written  in  the  following 
self-explanatory  shorthand  notation: 


(45>,<S4,=,S  £  *•••■  <"> 

»*• 


In  a  «Hrwii«r  one  three  more  useful  if  near  combinations: 


(45)  -  (54)  *  2  T[ 

»rt 

•  44 

•  *« 


(55)  ♦  (5«)  *  2  £ 

•  II 

a ad 


(55)  -  (56)  *  2  2 


« 


Now  baser  combinations  are  termed  between  Sqs.  (51)  and  (59)  or  between  Eqs.  (96)  sad 
(60)  with  the  purpose  of  reinforcing  or  canceling  terms  with  even  or  odd  »'i.  One  obtains 
tour  liasnr  comblsattoos  of  interest: 

(S»)  *  (»)  *  (45)  ♦  (54)  »  (55)  •  (56) 

•L 

HI 
#*•# 


i 
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(58)  f  (60)  s  (45)  -  (54)  4  (55)  -  (56) 


L  L 

•  II  •  1  1 
odd  even 

•  *  *  n’t 


(57)  -  (59)  =  (45)  4  (  54)  -  (55)  -  (56) 


L  L 

•  I  I  til 

even  odd 


c  ...  . 


(58)  -  (60)  *  (45)  -  (54)  -  (55)  4  (  56) 


L  L 

•  II  til 
odd  odd 

•  * l  n’t 


c  ...  . 


These  linear  combinations  are  expressions  exactly  of  the  form  needed  to  express  the 
amplitudes  in  Bq.  (53)  in  terms  of  eia. 

From  Bq.  (61)  one  obtains  the  amplitude  of  target  1  (strongest  target): 

Target  1  amplitude  *  0.33 

[/•(*»»  *,•>  ,  .'•<*,•**  *,.>  ,  t  («5) 

From  Bq.  (63)  one  obtains  the  amplitudes  of  the  second -strongest  true  target: 

Target  2  amplitude  -  0.35 

*i*>  .  *,*>  ,  ,>•<*„  *„••>  .  vi«tdl9**-dM*e)J  (M) 

The  smallest  true  target  amplitude  is  obtained  from  Bq.  (63): 

Target  3  amplitude  »  0.35 

Jrwtd,*-d„»  ,  ,«<*»,♦■*>„>  .  (57) 

The  false  target  amplitude  follows  from  Sq  (64): 

False  target  ampWade  »  0.25 

.  ,*•<*,*  *»•**>  ,  •  *„•*)]  (W) 

Iquatloes  (66)  through  (66)  are  complex  notations  deecrMag  both  the  ampHtuda  change, 
Le.,  capture  effect,  and  the  phase  change  due  to  limiting.  A  problem  exists  hew  to  pre- 
seat  the  smpltude  dlstrlbations  which  depeod  on  two  variables  (*„  and  *„)  and  on  two 
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parameters  (al  and  a3).  ft  was  decided  to  write  a  program  for  the  NRL  C DC -3800 
computer,  using  separate  real  and  imaginary  expressions  for  exp  ja  (»,.«,)  derived 
from  Ey.  (44): 


Re  e 


v#i> . 


1  ♦  a,  cos  0 j  ♦  cos  9j 


/<fe,  cos  0,  ♦  cos  )J  +  ( a,  sin  sin  )J 


(89) 


and 


7. 


*j  sin  S,  ♦  a j  sin  Sa 

^(1  ♦  cos  Sj  ♦  *j  cos  Sj  )*  ♦  (a,  sin  $t  *  *a  sin  )* 


(TO) 


The  product  obtained  from  the  computer  is  a  punched  paper  tape  which  contains  the  data 
for  a  pseudo-three -dimensional  plot  to  be  drawn  by  the  Gerber  plotting  machine  at  NRL. 
Thus  any  handling  of  numerical  data  by  a  draftsman  is  circumvented. 

This  report  shows  the  resulting  absolute  value  amplitude  distributions  for  three 
cases:  Figs.  11  are  for  three  equally  strong  signals  at  the  limiter  iiqwt;  Figs.  12  are 
for  two  equally  strong  signals  at  the  limiter  input,  with  the  third  signal  10  dB  smaller 
than  either  one  of  the  stronger  signals;  and  Figs.  13  are  for  one  strong  signal  and  two 
smaller  signals,  each  of  which  is  at  a  -10-dB  level  with  respect  to  the  strongest  signal. 
The  strongest  false  target  is  observed  if  the  three  true  target  signals  are  equally  strong 
and  in  phase  or  180  degrees  out  of  phase.  In  this  case  each  of  the  true  target  waveforms 
and  the  false  target  will  appear  with  an  amplitude  equal  to  1/2  of  the  maximum  possible 
amplitude.  With  Unite  length  codes  the  false  target  may  actually  be  slightly  stronger 
than  any  one  of  the  true  targets.  This  effect  is  due  to  the  Interference  of  range  aide 
lobes,  as  will  be  explained  shortly. 

U  one  hss  two  equally  strong  signals  and  one  smaller  signal  at  the  limiter  input, 
then  a  false  target  is  generated  of  approximately  the  same  aiae  as  the  captured  smaller 
signal,  aa  can  be  seen  from  comparing  the  surface  in  Fig.  12(d)  with  the  surface  in  Fig. 
12(c).  The  two  stronger  signals  suffer  capture  to  an  extent  depending  on  their  relative 
phiuiif^,  and  it  should  be  noted  (Figs.  12(a)  and  12(b))  that  each  maximum  of  one  signal 
coincides  with  a  minimum  of  the  other  signal  component. 

ft  has  been  discussed  that  there  is  no  false  target  generation  effect  if  there  are  only 
'wo  signals  at  the  limiter  input,  ft  should  therefore  be  no  surprise  that  the  iafoe  target 
«U*  la  much  reduced  If  two  of  the  three  true  target  signals  are  relatively  small  (Pig. 

3(d)).  The  discussion  of  capture  and  effects  of  false  target  generation  is  much  aimpli- 
ticd  if  all  input  signals  are  equally  strong  and  if  their  mutual  RF  phases  are  equal  to  0 
or  180  degrees.  Then  one  caa  read  from  Kq.  (44)  or  from  Bqs.  (80)  and  CTO)  that 
*»!>*«(«  .•,)  *  i  if  at  leest  one  of  the  ergur  mts  is  equal  to  0  and  that  <~*p  ;•<«.«)  *  -i. 
To  obtain  the  captured  target  and  the  false  target  amplitudes,  one  has  to  combine  the  ex¬ 
ponential  functions  according  to  Sqs.  (88)  through  (88).  One  obtains  the  following: 


Casel:  All  true  targets  having  the  same  polarity,  l.e.,  «. 


i« 


o: 


target  1  amplitude  >  0.28  ((1)  *  (1)  ♦  (1)  ♦  (-1)]  •  0.5 

target  2  amplitude  •  0.25  ((1)  -  (1)  ♦  (1)  -  (-i)J  *  0.5 

target  3  smpttude  •  0.25  ((1)  ♦  (I)  -  (1)  -  (-1)]  •  0.5 

•  0.25(d)  -  d)  -  d)  ♦  <-l>]  * 
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Fig.  11(a)  -  True  target  1  amplitude  for  the  case  of  three  equally  strong  expanded 
signals  using  zero-pi  shift  register  phase  coding  (for  machine  convenience  the 
parameters  <A10,  <i30,  a,,  and  a2  are  labeled  as  shown) 


Fig.  11(b)  -  True  target  2  amplitude  for  the  case  of  three  equally  strong  expanded 
signals  using  zero-pi  shift  register  phase  coding 
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TARGE <  NO  3 
AMPLITUDE 


u 


0 

0 


PH!  10  IDEG) 


PH!  20  (0EG1 


PARAMETER  CHOICE 
Rl= 1.00  R2= 1 o 00 


Fig.  11(c)  -  True  target  3  amplitude  for  the  case  of  three  equally  strong  expanded 
■Ignal ■  using  zero-pi  shift  register  phase  codirg 


Fig.  11(d)  -  False  target  amplitude  for  the  case  of  three  equally  strong 
expanded  signals  using  sero-pl  shift  register  phase  coding 
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Fig.  12<a)  -  True  target  1  amplitude  for  the  case  of  two  equally  strong  targets 
and  one  smaller  true  target  using  zero-pi  shift  register  phase  coding 


Fig.  12(b)  -  True  target  2  amplitude  for  the  oaae  of  two  equally  strong  targets 
and  one  smaller  true  target  using  zero-pi  shift  register  phase  coding 
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Fig.  12(c)  -  True  target  3  amplitude  for  the  case  of  two  equally  strong  targets 
and  one  smaller  true  target  using  zero-pi  shift  register  phase  coding 


Fig.  12(d)  -  False  target  amplitude  for  the  case  of  two  equally  strong  targets 
and  one  smaller  true  target  using  zero-pi  shift  register  phase  coding 
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Fig.  13(a)  -  True  target  1  amplitude  for  the  case  of  one  strong  target  and  two 
equal  Intensity  smaller  true  targets  using  zero-pi  shift  register  phase  coding 


TARGET  NO  2 
AMPLITUDE 


A 


Fig.  13(b)  -  True  target  2  amplitude  for  the  case  of  one  strong  target  and  two 
equal  intensity  smaller  true  targets  using  zero-pi  shift  register  phase  coding 
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TARGET  NO  3 
AMPLITUDE 
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PHI  20  IDEGI 


180 

PARAMETER  CHOICE 
A1  =  Od  32  A2=  0»  32 


Fig.  13(c)  -  True  target  3  amplitude  for  the  case  of  one  strong  target  and  two 
equal  intensity  smaller  true  targets  using  zero-pi  shift  register  phase  coding 
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Fig.  13(d)  -  False  target  amplitude  for  the  case  of  (me  strong  target  and  two 
equal  intenalty  smaller  true  targets  uaing  zero-pi  shift  register  phase  coding 
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After  hard  limiting  there  will  be  four  signals,  each  of  them  having  1/2  the  maximum 
possible  amplitude,  three  of  the  signals  corresponding  in  time  shift  and  polarity  to  the 
true  input  signals  and  the  fourth  signal  corresponding  to  a  false  target  with  a  polarity 
opposite  to  the  true  target  polarity. 

Case  2:  One  of  the  true  targets  having  phase  opposition  with  the  remaining  pair  of 
true  targets,  e.g.,  <t>i0  -  o  and  <620  =  v. 

target  1  amplitude  =  0.25  [(1)  +  (-1)  +  (1)  +  (1)J  =  0.5 

target  2  amplitude  =  0.25  [(1)  -  (-1)  +  (1)  -  (1)]  =  0.5 

target  3  amplitude  =  0.25  [(1)  +  (-1)  -  (1)  -  (1)]  =  -0.5 

false  target  amplitude  =  0.25  [(1)  -  (-1)  -  (1)  +  (1)]  =  0.5. 

The  limiter  output  contains  again  three  signals  of  the  same  polarity  and  one  signal  of  the 
opposite  polarity.  One  of  the  three  same -polarity  signals  corresponds  to  a  false  target, 
and  the  opposite  polarity  signal  corresponds  to  one  of  the  true  targets  in  this  case.  It 
should  be  obvious  that  the  polarity  of  the  compressed  signals  can  consequently  not  be 
evaluated  to  permit  discrimination  between  true  and  false  targets. 

Consider  the  output  signal  obtained  from  a  matched  filter  or  a  pulse  compression 
device.  In  the  single-target  case  one  will  obtain  an  impulse-like  point  target  response, 
consisting  of  a  central  spike  of  height  N,  where  n  is  equal  to  the  length  of  the  code,  sur¬ 
rounded  by  uniform -range  side  lobes  at  the  -1  level.  This  type  of  autocorrelation  func¬ 
tion  is  depicted  in  Fig.  14(a).  If  there  are  three  input  signals  of  equal  intensity,  corre¬ 
sponding  to  one  of  the  two  previously  considered  cases,  there  will  be  three  correlation 
spikes  of  one  polarity,  corresponding  to  three  true  targets  in  case  1  or  to  two  true  tar¬ 
gets  and  one  false  target  in  case  2,  of  a  magnitude  it/2  -  1,  and  one  correlation  spike 
with  the  opposite  polarity,  corresponding  to  a  false  target  or  a  true  target  and  with  a 
magnitude  it/ 2  >  l.  The  addition  or  subtraction  of  1  from  N/2  originates  from  the  con¬ 
structive  or  destructive  interference  of  a  uniform -range-side-lobe  level  with  the  corre¬ 
lation  peaks.  Hence  it  may  happen  that  the  false  target  is  larger  than  any  one  of  the  true 
target  responses. 

The  existence  of  a  spurious  target  response  has  been  confirmed  by  a  computer  sim¬ 
ulation.  Figure  14(a)  shows  the  autocorrelation  function  of  a  maximum -length  linear 
shift  register  code  with  l  ■  63  bits.  This  function  is  descriptive  of  the  radar  response 
to  a  single  point  target.  Figure  15(a)  is  a  plot  of  the  processor  response  due  to  three 
equally  strong  equi-phase  true  targets.  The  response  is  obtained  through  addition  of 
three  mutually  delayed  binary-coded  video  signals,  hard  limiting,  and  crosscorrelation 
with  the  original  code.  Hie  resulting  crosscorrelation  function  shows  a  spurious  re¬ 
sponse  exactly  as  predicted  by  the  theory.  The  preservation  of  a  uniform -range -side  - 
lobe  level  and  the  predicted  capture  effect  are  confirmed.  The  uniform-range-side-lobe 
level  is  preserved  since  the  crosscorrelation  function  depicted  in  Fig.  15(a)  is  the  same 
as  the  sum  of  four  autocorrelation  functions,  each  of  them  having  uniform  side-lobe  lev¬ 
els  as  shown  in  Fig.  14(a),  weighted  with  the  factors  +1/2  or  -1/2,  and  mutually  delayed 
by  the  amounts  Tlt  r„  and  r. 

Figure  16(a)  when  compared  with  Fig.  15(a)  exhibits  that  the  false  target  may 
undergo  a  very  erratic  motion  if  the  true  target  geometry  changes  but  very  slightly. 

This  erratic  motion  effect  may  be  used  in  a  video  integration  process  to  enhance  the  ra¬ 
tio  of  true  targets  to  false  targets. 
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Fig.  14(a)  -  Simulated  radar 
response  to  a  single  target, 
for  a  radar  using  a  63-bit 
repetitive  code 


Fig.  14(b)  -  Simulated  radar  response  to  a  single 
target,  for  a  radar  using  a  63-bit  truncated  code 


Fig.  15(a)  -  Simulated  radar 
response  to  three  point  targets, 
for  a  radar  using  a  63-bit  re¬ 
petitive  code 
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Fig.  15(b)  -  Simulated  radar  response  to  three  point 
targets,  for  a  radar  using  a  63-bit  truncated  oode 
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Fig.  16(a)  -  Simulated  radar 
response  to  three  point  targets, 
of  a  slightly  varied  geometry, 
for  a  radar  using  a  63-bit  re¬ 
petitive  code 


Fig.  16(b)  -  Simulated  radar  response  to  three  paint 
targets,  of  a  slightly  varied  geometry,  for  a  radar 
using  a  63-bit  truncated  code 


Figures  14(b),  15(b),  and  16(b)  show  the  poise  compression  waveforms  obtained  if 
the  code  sequences  are  truncated.  It  is  known  that  truncation  increases  the  average 
range-side-lobe  level.  The  spurious  targets  may  be  observed  at  the  same  locations 
where  they  would  be  if  the  code  sequence  were  repeated.  The  false  target  response  may 
be  of  reduced  magnitude  and  depends  on  how  much  the  expanded  radar  signals  overlap. 
Since  the  range-side-lobe  structure  of  truncated  sequences  Is  very  complicated  and  ir¬ 
regular,  an  exact  theory  to  quantitatively  explain  capture  and  false  target  generation  ef¬ 
fects  should  not  be  expected. 


CONCLUSIONS 

ft  is  shown  in  this  report  that  hard  limiting  of  superimposed  expanded  radar  returns 
may  cause  a  false  target  effect  in  addition  to  the  well-known  capture  effect.  The  false 
target  effect  has  been  known  for  the  case  of  a  pair  of  linear  Fll  radar  returns,  ft  is 
shown  that  a  zero-pi-phase-coded  radar  signal  modulated  by  a  maximum-length  shift 
register  sequence  may  be  subject  to  a  false  target  effect,  ft  is  interesting  to  note  that 
with  two  targets  in  the  case  of  sero-pi  phase  coding  there  is  capture  only  and  no  false 
target  effect;  at  least  three  true  targets  must  be  taken  into  account  to  show  the  existence 
of  a  false  target.  This  study  may  serve  as  an  example  to  show  the  erroneous  conclusions 
one  may  reach,  in  nonlinear  systems  theory,  by  deriving  rules  about  multiple  target  re¬ 
sponses  from  the  study  of  a  two-target  model. 

To  avoid  false  target  effects  <mw  must  avoid  the  limiting  of  expanded  radar  returns 
having  a  high  signal -to-nolse  ratio. 
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Appendix 

NUMERICAL  DETERMINATION  OF  THE  FOURIER  COEFFICIENTS 


The  coefficients  c„  are  the  Fourier  coefficients  of  the  complex  function  «  'p  ;  u " ) 
as  can  be  seen  from  Eq.  (8).  The  List  Fourier  transform  is  a  convenient  technique  for 
the  given  task,  in  particular  since  it  is  "naturally”  adapted  to  the  handling  of  complex 
variables.  Therefore  a  program  ha*  been  written  in  the  XTRAN  language  for  a  time¬ 
sharing  computer  terminal  facility  at  NRL,  using  an  FFT  subroutine  provided  by  Com- 
share.  A  program  listing  is  reproduced  at  the  end  of  this  appendix. 

The  program  goes  through  the  following  logical  steps: 

1.  Read  *  «  number  of  samples;  this  must  be  a  power  of  2. 

2.  Read  the  small -to-largc  signal  ratio. 

3.  Compute  a  period  of  the  function  .xp  )i<»  j  at  the  *  sampling  points 


using  separate  numerical  expressions  for  the  real  and  .'or  the  imaginary  parts  of  the 
function 

„>•<“>  I  <*»’*  U  l  «if  <)|  (  %•  a  0 

I  •  .i  ro*  f>  n  ttn  n 

/r-  .I*  j.t  «*«-.  >  yr  *  »»* >» 

4.  Calculate  the  Fourier  coefficients  of  this  function  by  means  of  an  FFT  subroutine. 

5.  Print  out  all  r„  values  of  (Merest  (the  printout  restricted  to  the  range 
«  •  -13,  -14,  .. . ,  14,  15  and  c„  -80  dB). 

6.  Go  to  step  1  to  reiidttate  the  program  with  a  new  data  set.  (Subsequent  to  step  5 
the  program  returns  to  its  starting  point  so  that  a  new  data  set  can  be  read  and  the  cor¬ 
responding  c„  values  can  be  determined.  The  program  recycles  until  its  execution  is 
interrupted  by  an  "escape*'  command  to  the  computer.) 

To  ensure  that  a  sufficiently  large  number  of  sampling  points  was  used,  their  num¬ 
ber  was  increased  by  factors  of  2  from  one  iteration  to  the  next  until  the  calculated  co¬ 
efficients  rem.:  ned  the  same  to  within  0.01  dB.  To  exhibit  the  convergence  of  the 
computed  <'»  values  toward  a  limit  when  the  number  of  samples  is  increased,  the  se¬ 
quence  of  printouts  for  N  ■  *  -10  dB  and  for  »  «  8,  is,  32,  and  128  samples  Is 
reproduced  within  Tabic  At. 
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Table  A1 

Values  Printed  Out  for  20  log  a  •  -10  dB 
and  for  *»8,  16,  32,  and  128  Samples 


Order 

C„  (dB) 

Order 

Cn 

Cn  (dB) 

Number  of  Points  it  •  8 

Number  of  Points 

*  «  32 

-4 

0.0022 

-53.14 

-4 

0.0026 

-51.67 

-3 

-0.0094 

-40.58 

-3 

-0.0094 

-40.50 

-2 

0.0359 

-28.90 

-2 

0.0359 

-28.90 

-1 

-0.1521 

-16.36 

-1 

-0.1521 

-16.36 

0 

0.9745 

-0.22 

0 

0.9745 

-0.22 

1 

0.1601 

-15.91 

1 

0.1602 

-15.91 

2 

-0.0126 

-37.98 

2 

-0.0128 

-37.84 

3 

0.0013 

-57.72 

3 

0.0020 

-53.80 

Number  of  Points  it  »  16 

Number  of  Points  *  »  128 

-4 

0.0026 

-51.67 

-4 

0.0026 

-51.67 

-3 

-0.0094 

•40.50 

-3 

-0.0094 

-40.50 

-2 

0.0359 

-28.90 

-2 

0.0359 

-28.90 

-1 

•0.1521 

-16.36 

-1 

-0.1521 

-16.36 

0 

0.9745 

-0.22 

Q 

0.9745 

-0.22 

1 

0.1602 

-15.91 

1 

0.1602 

-15.91 

2 

-0.0128 

-37.84 

2 

-0.0128 

-37.84 

3 

0.0020 

-53.80 

3 

0.0020 

-53.80 

In  UUs  case  Id  cfeda  points  are  sufficient.  The  printouts  remain  the  same  for  a 
larger  number  of  data  points,  as  shown  for  *  -  32  and  128.  The  use  of  only  eight  data 
points  would  entail  some  error,  especially  as  for  as  the  relatively  smaller  else  beat 
products  are  concerned. 

The  largest  number  of  samples  was  required  in  the  vicinity  of  a  1 ,  i.e.,  if  *  ) 

was  close  to  a  linear  sawtooth  function.  In  this  caae  a  1024 -point  transform  was  selected. 
For  smaller  a  values  s  much  smaller  number  of  sampling  points  sufficed. 

The  validity  of  the  program  was  checked  by  the  following  criteria: 

1.  The  imaginary  part  of  the  Fourier  coefficients  turned  out  to  be  0  within  the  con* 
fines  of  roowk.fi  errors,  as  required  by  the  theory  presented  In  the  section  "General 
Theory  of  the  Hxrd  Uniting  of  Two  Input  Signals," 

2.  For  •>  1  (s,  s,  ■  0  dB)  the  coefficients  obtained  by  the  computer  program 

agree  vita  those  calculated  from  the  spec lal-case -equation  (IS)  to  within  about  0.01  dB. 

3.  For  small  -•  values  the  6.02-dz>  loss  as  suggested  by  Sq.  (10)  has  been  confirmed. 

4.  For  all  values  of  the  calculated  c„  values  are  In  excellent  agreement  with 
experimentally  obtained  values.  The  experiment  consists  of  the  spectral  analysis  of  the 
hard* limited  summation  of  two  sinusoids,  as  described  in  the  section  "Limiting  of  a  Pair 
of  Constant-Frequency  CW  Sinusoids." 

The  computer  results  are  presented  in  two  formats: 
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1.  A  set  of  curves  for  cr.  (  n)  with  n  values  between  >4  and  5  and  with  <•<  values  cor¬ 
responding  to  the  range  from  0  to  -30  dB,  in  Fig.  6. 

2.  Line  spectra  including  a  numerical  printout  for  the  c„  values  in  Fig.  7. 

Both  formats  have  their  specific  merits  and  drawbacks.  The  advantage  of  format  1 
is  to  indicate  clearly  how  the  c„  values  vary  as  a  function  of  The  drawing  would  be 
too  crowded  if  one  would  attempt  to  display  in  it  all  higher  order  beat  products,  however. 
For  this  reason  the  presentation  in  format  2  was  added.  Figure  7  shows  all  beat  prod¬ 
ucts  whose  order  is  between  -14  and  15  if  they  are  stronger  than  -50  dB  with  respect  to 
the  total  signal  energy. 


NRL  REPORT  7001 


XTRAN  Program  for  the  Determination  of  cn(a>; 

D’N  ARRAY(2, 1024, 1, I ), 1 ( 3) , M(*s6> * SAkkAY(2S6) 

PI  =  3.I41S926535 
NPG=B 

*  NIP  =  TOTAL  NUMBER  0F  UATA  P0IN1 S,  MAY  BE  Sl'BDl VIDEO  INTO 
NC>  GROUPS  WITH  NPG  1)4 f A  POINTS  PER  GROUP 
I  WRITE  (1,11) 

READ  10,32)  NT P 
NG  =  NTP/NPG 

♦EVALUATE  THE  FORMULA  FOR  ALPHA 
WRITE  (1,13) 

READCO, 1 7)  AIA 

A  =  EXP(AL0G(1O« )/20«*(“AIA)) 

SI  «  SIN(2.*PI/NTP) 

CO  =  C0S(2.*PI/NTP) 

DU  200  J=1,NG 
S  =  SIN(g.*PI/NG*<J-l  )  ) 

C=  C0S(2.*PI/NG*(J-I > ) 

00  100  LL=I,NPG 
IF  (LL<EO. 1 )  GOTO  SO 
STEMP  *  SI  *c+ca  *S 
C  =  C0  *C-  SI  *3 
S  =  STEMP 

SO  DEN  =  SORT ( I  A*(A+2.*C>) 

ARRAY( 1 ,NPG*( J“ I )+LL, 1,1)  *  (l.+A*C)/DEN 
100  ARRAY(g,NPG*(J-l >*LL,1, 1 )=A*S/DEN 
300  C’E 
WRITE  (1,14) 

ANfPsNTP 

1(1 )=  AL0G(ANTP)/AL0G(2«  >  ♦ 

FOR  L=2,3i  I(L)=0 

CALL  F A SFT ( ARRA Y, I , M, SARRAY, - 1 , I FL ) 

WRITE  (1,1S) 

00  SOO  N=MAXO (NTP/2+ 1 ,NTP- 1 S),NTP 
X=ARRAY(|,N, 1,1) 

IF  (AB5(X).LT..001 )  G0T0SOO 

WRITE (1,16)  N-1-NTP,X, 10. AAL0G10(X*X) 

500  C’E 

D0  600  N=l,MIN0(NTP/2, 15) 

X=  ARRAYd  ,N,  1,1  ) 

IF(ABS(X).LT. .001 )  GOT06OO 
WRITE! 1, 16)  N-1,X, 1O.*AL0G1O(X*X> 

600  C’E 
WRITE  (1,14) 

G'fl  1 

11  F ’T(3/,$SPECIFY  TOTAL  NUMBER  OF  POINTS  I  S,Z> 

12  F’T  (14) 

13  F’T  (3X,i  A  [MINUS  D8)  1  $,Z> 

14  F’T  (3/) 

is  F’T  (/,SORDER  CN  CN  IN  DBS, 3/) 

16  F*T(X, I3,4X,F10.4,5X,F9.2) 

17  F’T(F5.2) 

END 
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13.  ABSTRACT 


Hard  limiting  before  pulse  compression  or  correlation  processing  is  a  common  approach 
to  the  CFAR  (constant  false  alarm  rate)  problem,  and  it  offers  a  good  and  simple  solution  in  a 
single-target  or  scarce-target  environment.  With  the  advent  of  radars  with  a  large  time- 
bandwidth  product  the  possibility  arises  that  expanded  radar  returns  due  to  multiple  targets 
of  interest  may  overlap  very  largely  or  entirely  but  still  may  be  sufficiently  separated  to  be 
resolved  after  receiver  processing.  In  this  case  the  compressed  pulses  cannot  attain  full 
amplitude  at  the  processor  output  even  if  the  signal-to- noise  ratio  at  the  input  is  very  high; 
this  phenomenon  is  known  as  capture  and  small  signal  suppression.  The  purpose  of  this  re¬ 
port  is  to  exhibit  that,  in  addition  to  compressed  target  responses  of  reduced  magnitudes, 
false  targets  may  be  generated  with  apparent  amplitudes  of  the  same  order  or  exceeding  those 
of  legitimate  targets.  Spurious  target  generation  in  the  case  of  chirp  radar  has  been  known 
for  some  time.  The  theory  has  been  extended  to  maximum-length  linear  shift- register  codes 
which  are  used  as  modulation  functions  of  pulse-compression  and  phase-coded  CW  radars.  It 
is  found  that  a  single  pair  of  radar  returns  coded  in  this  manner  is  subject  to  capture  only  and 
not  to  false  target  generation.  Surprisingly,  however,  the  addition  of  a  third  expanded  signal 
produces  a  spurious  response.  This  generation  of  a  false  target  should  be  taken  in  account 
when  the  dynamic  range  of  future  phase-coded  radars  using  linear  shift- register  codes  is 
specified,  in  particular  if  the  radar  is  designed  for  automatic  track  and  raid-size  determina¬ 
tion. 
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